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Abstract. A 2.5-dimensional magnetohydrodynamics simulation analysis of the energy release for three different reconnec- 
tion regimes is presented. The system under investigation consists in a current-sheet located in a medium with a strong density 
variation along the current layer: such system is modeled as it were located in the high chromosphere/low solar corona as 
in the case of pre- flare and coronal mass ejection (CME) configurations or in the aftermath of such explosive phenomena. 
By triggering different magnetic-reconnection dynamics, that is from a laminar slow evolution to a spontaneous non-steady 
turbulent reconnection [1, 2, 3], we observe a rather different efficiency and temporal behavior with regard to the energy fluxes 
associated with each of these reconnection-driven evolutions. These discrepancies are fundamental key-properties to create 
realistic models of the triggering mechanisms and initial evolution of all those phenomena requiring fast (and high power) 
magnetic reconnection events within the solar environment. 
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INTRODUCTION 

Magnetic reconnection is essential to model partially or 
completely most of the phenomena occurring within the 
solar atmosphere, being the sun's magnetic field the huge 
reservoir of energy driving the dynamical evolution of 
such structures. A detailed and realistic model of magne- 
tohydrodynamics (MHD) reconnection is an unescapable 
step in order to achieve a deep understanding of the 
highly dynamic solar atmosphere as well as of the in- 
teraction of solar plasma structures with the interplane- 
tary medium. In fact, the MHD approach is adequate to 
describe the initial evolution and the following nonlin- 
ear stages of such phenomena, even though the detailed 
physics of their initiation depends primarily on kinetic 
effects. 

In general, the starting point is either the well-known 
Sweet-Parker's model of a slow and laminar reconnec- 
tion [4, 5] in a harris or force-free current-sheet, or the 
fast process described by Petschek's model [6] and usu- 
ally obtained by the introduction in the MHD system of 
a localized enhancement of the plasma resistivity [7, 8]. 
Yet, a new high-power and non-steady turbulent recon- 
nection process is possible in MHD without invoking 
the presence of any ad-hoc term resembling kinetic ef- 
fects in MHD equations. This has been recently shown 
by Lapenta [1] and applied to the solar environment 
by Bettarini and Lapenta [2], Such mechanism includes 
a spontaneous transition from a slow to a fast reconnec- 
tion with time-scales of the order of Alfven times and no 
factitious anomalous effects are considered. During its 



evolution, current-sheets develop a chaotic structure of 
multiple small-scale interacting reconnection sites deter- 
mining a macroscopic turbulent state. 

As in Bettarini and Lapenta [2], the system under in- 
vestigation in the present work is a current-sheet embed- 
ded in the high chromosphere/low solar corona region. 
It consists in a current-sheet determined by a force-free 
magnetic field in the presence of a strong density vari- 
ation along the current layer provided by a proper step 
function. We vary the global resistivity of our numeri- 
cal domain (by means of the Lundquist number) such 
that different dynamical regimes of the current-sheet are 
set, from a slow and diffusive reconnection to a fast and 
turbulent process and an in-between evolution. Here, we 
present an energy-release analysis of such three differ- 
ent processes in order to (a) provide useful hints of the 
best-promising mechanism to be considered for more 
and more realistic models of solar (explosive) phenom- 
ena, and (b) to determine the key-features of events that 
could be observed by current and near-future solar and 
space missions. 

EXPERIMENT SETTINGS 

We solve numerically the viscous-resistive compress- 
ible 2.5D one-fluid MHD equations parametrized by the 
global Lundquist and kinetic Reynolds number, S and R v 
respectively. We define a vertical or stream-wise direc- 
tion (z) oriented away from the sun, wherein open bound- 
ary conditions are set at z — and L z , and a horizon- 
tal or cross-stream direction (x) where we set reflecting 



Energy release into the solar atmosphere from magnetic reconnection 



September 20, 2009 



1 



Table 1. Summary of the simulations 
performed for the present work and the 
related parameters. 

Run A Run B Run C 

S 10 2 10 3 10 4 

j8 = 0.2, R v = 10 4 ,y=5/3 
xe [-10.,10.],z€ [0.,80.] 
dx = 0.1, dz = 0.25 



boundary conditions. We consider an ideal equation of 
state with poly tropic index y such that p = p (y — 1)1 
where p, p and / are respectively the kinetic pressure, the 
density and the enthalpy of the plasma. We use the La- 
grangian code, FLIP3D-MHD [9], and we consider 600 
(in x) x 960 (in z) Lagrangian markers arrayed initially 
in a 3 x 3 uniform formation in each of the 200 x 320 
cells of our numerical grid. As already pointed out in 
the introduction, our system consists in a current-sheet 
configuration determined by a force-free magnetic field 
with B z (x) — tanhx, B x (x) = and the out-of-plane com- 
ponent of the magnetic field B y (x) = — sechx as guide 
field. With this configuration, we obtain high-resolution 
simulations where the current-sheet width is resolved by 
30 Lagrangian markers. The density is modeled as a step 
function where the ratio between the density in the region 
< z < L z /% (modeling the high chromosphere) and that 
one in L z /8 < < L z (the low corona) is about 10 5 as 
considered in Yokoyama and Shibata [8], Bettarini and 
Lapenta [2]. We consider an initial GEM-like perturba- 
tion [10] but exponentially localized at the center of the 
numerical box (i.e. within the solar corona) with a per- 
turbation amplitude of e = 0.5. 



the global resistivity) is shown. This quantity is mea- 
sured by computing the maximum and minimum val- 
ues of the out-of-plane magnetic potential, A y at the mid 
axis (x — 0): their difference provides the amount of flux 
in the closed field lines, caused by reconnection as de- 
scribed in Biskamp [7]. A low value of S, S= 10 2 results 
in a laminar and slow dynamics of the current-sheet. Af- 
ter an initial transient time, a nearly-constant reconnected 
flux is observed (Run A, dashed line). By setting S to 
10 3 , the system presents again an initial nearly-constant 
increase in the reconnected flux resembling a slow re- 
connecting evolution, although it is more and more evi- 
dent a smooth transition to a faster evolution as the time 
goes by (Run B, dotted line). By increasing S further, 
S = 10 4 , smaller scales can be reached and the high 
resolution allows us to observe a completely different 
regime [1,2]. Because of a tearing instability of the ini- 
tial Sweet-Parker configuration [1 1, 3], the system expe- 
riences an abrupt switch from an initial slow and laminar 
reconnection regime to a fast reconnection process lead- 
ing to a macroscopic turbulent configuration of the initial 
system. The different behavior of the reconnection pro- 
cess in the three simulations is pointed out also in the 
evolution of the reconnection rate as a function of time. 
This is shown in the bottom panel of Fig. 1 . Different 
rates but the same general behavior are shown during the 
first few decades of time steps where the system is driven 
by a slow reconnection. At about t = 40 (measured in 
terms of alfven times), it is evident the rapid increase of 
the reconnection rate in Run C. With regard to this case, 
the rate is measured both by computing the slope of the 
corresponding reconnected flux (solid line in Fig. 1) as 
for the other runs and also by means of the out-of-plane 
electric field variation through the reconnection diffusion 



Physical Parameters 

The MHD equations are normalized according to 
reasonable values for the high chromosphere/low solar 
corona: the number density p/m 10 9 cm~ 3 ; the length 
scale L w 1000 — 2000km, where L is the current-sheet 
width and L = L x /20 = L,/80; the Alfven velocity « 
400 — 600 kms -1 (approximately three times the sound 
speed). So, we have a time scale of = L/va ~ 1 -5 — 5 s. 
See Tab. 1 for more details on simulations' setting and 
parameters. 

RESULTS 

In the top panel of Fig. 1, the reconnected magnetic 
flux (in dimensionless units) as a function of time for 
three different values of the Lundquist number (so, of 
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Figure 1. Top panel: Measured reconnected flux as a function 
of time where the dashed line corresponds to Run A, dotted line 
to Run B and solid and dotted-dashed lines correspond both 
to Run C (see Tab. 1). Bottom panel: Reconnection rate as a 
function of time for the three simulations. 
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region (dotted-dashed line). With the latter method, the 
curve of the rate is more noisy, yet the two methods give 
two manifestly consistent results. Note that also the rate 
of Run B (dotted line) is becoming more and more con- 
sistent with a fast evolution as the time goes by, whereas 
the Run A's rate is slowly decreasing. 

As described in Bettarini and Lapenta [2], during such 
three different regimes, the slow reconnection evolution, 
the spontaneous and abrupt variation from a slow to a fast 
regime and the intermediate dynamics, the current-layer 
develops asymmetrically because of the strong density 
variation resembling the transition from the high chro- 
mosphere and the low solar corona. Moreover, character- 
istic features in magnetic field topology of the aftermath 
of solar explosive phenomena are observed. In order to 
understand the importance of such reconnection mech- 
anisms within the solar atmosphere, a fundamental step 
is to analyze such processes from the point of view of 
their efficiency in releasing energy and particles in the 
upper atmosphere. In Fig. 2, the flux of the Poynting 
vector, the enthalpy flux, the kinetic flux, and the mass 
flux through the upper boundary of the numerical box 
(z = L z ) are shown as functions of time. In the figure 
also the normalization values for all the shown quanti- 
ties are reported. As the perturbation is approximately in 
the middle of our numerical box (z ~ L z /2 = 40, that is 
in the solar corona) and there the alfven velocity is 1, we 
start to measure the effects of reconnection at the upper 
boundary from t ~ 40. Of course, we neglect the initial 
transient time caused by a perturbation that is not exactly 
an eigenstate of the system. 

In the Run A, the current-sheet evolves accordingly 
a diffusive reconnection process that determines a pile- 
up of the enthalpy flux of the system (top-left panel of 



Fig. 2) as well as an increase of the mass flux (top- 
right panel) to the upper layers of the corona because 
of an acceleration process as shown by the kinetic flux's 
profile (bottom-right panel of Fig. 2). The flux through 
the upper boundary of the stream-wise component of the 
Poynting vector is almost null throughout the simulation. 
As previously pointed-out, Run A's reconnection rate 
sensitively decreases in the last phase and after a plateau 
of a few time steps, we observe a corresponding decrease 
of all fluxes. 

The slow-to-fast reconnection evolution of the current- 
sheet observed in Run C passes through a chaotic state 
during which the formation of several smaller-scale 
current-sheets and reconnection sites interacting each 
other is observed [2]. So, it is the interaction of such 
structures continuously moving backward and forward 
that prevents the system to release a large amount of en- 
ergy and particles to the upper atmosphere. Yet, as the 
process goes on and the system reaches a macroscopic 
turbulent configuration, an impulsive event is observed 
and a release of mass as well as an increase of the temper- 
ature is rapidly determined in the upper corona. A differ- 
ent consideration has to be done with regard to Poynting 
vector's flux. In fact, it increases constantly and sensi- 
tively during all the fast regime, but it does not charac- 
terize the impulsive phase of the system. 

As already pointed out, an intermediate evolution is 
observed during the Run B. The system switches to 
a fast reconnection regime according to a slow rather 
than an abrupt process. From the energy- and mass-flux 
point of view, the profiles shown in 2 stay in-between 
the behavior described for the previous cases. An initial 
slow but constant increase of the fluxes is followed by 
a plateau of few time-steps that does not lead to a final 
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decrease of the fluxes as in the Run A but to a impulsive 
phase as in the Run C, though being the peak values 
sensitively lower. 

CONCLUSIONS 

Here, we observe the energetic and mass release of a 
force-free current-sheet evolving in a small region of the 
high/chromosphere/low solar corona. According to the 
new mechanism proposed by Lapenta [1] and applied 
by Bettarini and Lapenta [2] within the solar environ- 
ment, it is possible that this system follows a magnetic 
field-line reconnection dynamics and that spontaneously 
evolves from a low- and constant-rate phase to a fast 
and high power regime within the pure-resistive MHD- 
theory framework (and so without invoking the presence 
of any anomalous resistivity). Hence, it becomes funda- 
mental the understanding whether this process can pro- 
vide us with evolution-related observables that can im- 
prove our models of solar (explosive) phenomena and of 
MHD structures in general. 

By performing 2.5D MHD high-resolution simula- 
tions and by varying the characteristic global Lundquist 
number, we switch three different magnetic reconnec- 
tion regimes: a Sweet-Parker evolution, the new spon- 
taneously slow-to-fast reconnection and an intermediate 
dynamics. By means of the Sweet-Parker regime, we can 
obtain an quantitatively efficient but slow conversion of 
the initial magnetic energy available in the basic config- 
uration, but it does not provide some important features 
observed for example in solar phenomena like flares or 
CMEs, that is a spontaneous impulsive event that evolves 
on time-scales of tens/hundreds of seconds (see again 
Fig. 2 and the normalization values there reported) and 
releases large quantities of mass, enthalpy and high- 
velocity structures that eventually can power and trigger 
secondary events in the upper solar atmosphere. Even an 
intermediate evolution combining an initial laminar re- 
connection dynamics and evolving to a faster and faster 
regime is not able to recover such dynamical features, 
even if a less rapid and power impulsive process is ob- 
served. Here, we show that the spontaneous non-steady 
reconnection process can provide time-scales, reconnec- 
tion rates and energy- and mass-flux values that, on the 
one hand, are consistent with the theoretical framework 
of pure resistive-MHD and, on the other hand, they allow 
to use such mechanism as starting model of the complex 
triggering event and initial evolution of solar phenomena. 

Future studies will include a detailed analysis both 
from the dynamical and energetic point of view of the 
smaller-scale structures forming during the chaotic evo- 
lution of the current-sheet by means of adaptive-mesh re- 
finement numerical tools. Furthermore, at this stage the 
introduction of a more realistic model of the high chro- 



mosphere/low solar corona in a three-dimensional con- 
figuration is mandatory. 
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